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We study direct CP violation in the hadronic decay B^ -^ p^ir^, including 
Mh. the effect oi p — uj mixing. We find that the CP violating asymmetry is strongly 



dependent on the CKM matrix elements, especially the Wolfenstein parameter r/. 
For fixed Nc (the effective parameter associated with factorization), the CP violating 
asymmetry, a, has a maximum of order 30% — 50% when the invariant mass of the 
vr+vr^ pair is in the vicinity of the iv resonance. The sensitivity of the asymmetry, 
a, to Nc is small. Moreover, if Nc is constrained using the latest experimental 
branching ratios from the CLEO collaboration, we find that the sign of sin 5 is 
always positive. Thus, a measurement of direct CP violation in B -^ p^n would 



remove the mod(7r) ambiguity in arg 
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1 Introduction 

Even though CP violation has been known since 1964, we still do not know the source 
of CP violation clearly. In the Standard Model, a non-zero phase angle in the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix is responsible for CP violating phenomena. In the 
past few years, numerous theoretical studies have been conducted on CP violation in 
the B meson system |jl], 0]. However, we need a lot of data to check these approaches 
because there are many theoretical uncertainties - e.g. CKM matrix elements, hadronic 
matrix elements and nonfactorizable effects. The future aim would be to reduce all these 
uncertainties. 

Direct CP violating asymmetries in B decays occur through the interference of at least 
two amplitudes with different weak phase (j) and strong phase 6. In order to extract 
the weak phase (which is determined by the CKM matrix elements), one must know 
the strong phase 6 and this is usually not well determined. In addition, in order to 
have a large signal, we have to appeal to some phenomenological mechanism to obtain 
a large 6. The charge symmetry violating mixing between p° and u, can be extremely 
important in this regard. In particular, it can lead to a large CP violation in B decays 
such as -B^ -^ p°(u;)7r^ — > tt+tt^tt^, because the strong phase passes through 90° at the u 
resonance |^, ^, ^]. Recently, CLEO reported new data on 5 ^ pn. It is the aim of the 
present work to analyse direct CP violation in B^ -^ p^{uj)'K^ -^ tt+tt^tt^, including p~uo 
mixing, using the latest data from the CLEO collaboration to constrain the calculation. 
In order to extract the strong phase (5, we use the factorization approach, in which the 
hadronic matrix elements of operators are saturated by vacuum intermediate states. 



In this paper, we investigate five phenomenological models with different weak form factors 
and determine the CP violating asymmetry for B"^ -^ p^{uj)tt^ — > 7r^7r~7r^ in these 
models. We select models which are consistent with the CLEO data and determine the 
allowed range of Nc (0.98(0.94) < Nc < 2.01(1.95)). Then, we study the sign of sin 5 in 
the range of Nc allowed by experimental data in all these models. We discuss the model 
dependence of our results in detail. 

The remainder of this paper is organized as follows. In Section 2, we present the form 
of the effective Hamiltonian and the values of Wilson coefficients. In Section 3, we give 
the formalism for the CP violating asymmetry in B^ —>■ p^{uj)7i~^ -^ 7i~^7i~7i~^, for all the 
models which will be checked. We also show numerical results in this section (asymmetry, 
a, the value of sin (5). In Section 4, we calculate branching ratios for B^ -^ p^ir^ and 
B^ -^ p+7r^ and present numerical results over the range of N^ allowed by the CLEO 
data. In last section, we summarize our results and suggest further work. 

2 The effective Hamiltonian 

In order to calculate the direct CP violating aymmetry in hadronic decays, one can use 
the following effective weak Hamiltonian, based on the operator product expansion P, 
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q=d,s i=3 



^AB=i = ^[ E K.K;(cior + c^o^,) - Vt,v:^ y. ^.o.] + /^.c, (i) 



where Cj(i = 1, ■ ■ ■, 10) are the Wilson coefficients. They are calculable in renormalization 
group improved pertubation theory and are scale dependent. In the present case, we use 



their values at the renorinahzation scale /i ~ nib. The operators Oi have the following 
form, 

Oi = qal^,{l - 75)m/3M/37''(1 - 75)&a, O^ = g7^(l - 75)mm7'^(1 - 75)6, 

O5 = g7M(l - 75)& E ^V(l + 75)g', Oe = qal,{l - l,)h E ^'pl''^^ + Ts)?:,, 
9' 9' 

3 

07= 2^7m(1 -75)&Ee9'^V(l + 75)g', 
g' 
3 

Os = -ga7/.(l - 75)^/3 E eg'g^7''(l + 75)ga, 

q' 
3 

09= 2^7/.(l-75)&Ee9'^V(l-75)g', 

3 

Cio = 2^"7m(1 - 75)^/3 E eg'g^7^(l - 75)ga, (2) 

9' 

where a and (3 are color indices, and q' = u, d or s quarks. In Eq.(2), Oi and O2 are 
the tree level operators, O3 — Oq are QCD penguin operators, and O7 — Oiq arise from 
electroweak penguin diagrams. 

The Wilson coefficients, Cj, are known to the next-to-leading logarithmic order. At 
the scale /i = m;, = 5GeV, they take values the following values ||, |]: 

ci = -0.3125, C2 = 1.1502, 

C3 = 0.0174, C4 = -0.0373, 

C5 = 0.0104, C6 = -0.0459, 

C7 = -1.050 X 10-^ cg = 3.839 x 10-^ 

cg = -0.0101, cio = 1.959 X 10"^ (3) 



To be consistent, the matrix elements of the operators Oi should also be renormalized to 
the one-loop order. This results in the effective Wilson coefficients, c^, which satisfy the 
constraint, 

Q(m,)(0,(m,)) = c:(0,)*^'=^ (4) 

where (Oj) '^^^ are the matrix elements at the tree level, which will be evaluated in the 
factorization approach. From Eq.(4), the relations between c[ and Ci are p, ||, 



Ci — Ci, C2 — C2, 



C3 = C3 - Ps/3, C4 = C4 + Ps, 

C5 = C5 — Ps/3, Cg = Cg + Ps, 

Cy = C7 + Pe; Cg = Cg, 

Cg = Cg + Pe, c'lo = Cio, (5) 



where 



Ps = {as/Sn)c2{lO/9 + G{m„ /x, q^)), 

Pe = {aem/97T){3ci + C2)(10/9 + G{m„ fi, g2)), 



with 



G{mc, /i, g^) = 4 / dxx{x — l)ln 

JO 



1 m^ — a:(l — a;)g^ 



/^^ 



Here g is the typical momentum transfer of the gluon or photon in the penguin diagrams. 



Glrric, fj^jq"^) has the following explicit expression [10|, 
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(6) 



Based on simple arguments at the quark level, the value of g^ is chosen in the range 
0.3 < q^ /ml < 0.5 P, Q. From Eqs.(5,6) we can obtain numerical values for c[. 
When q'^/ml = 0.3, 



c[ = -0.3125, 4 = 1.1502, 



Co = 2.433 X 10^2 + 1.543 x lO^^i, c, = -5.808 x 10"^ - 4.628 x 10 



ci = 1.733 X 10^2 + 1.543 x lO^^z, c. = -6.668 x 10"^ - 4.628 x 10"^^ 



' "-6 



c'7 = -1.435 X 10"^ - 2.963 x lO'^i, Cg = 3.839 x 10" 



4 = -1.023 X 10^2 - 2.963 x lO^^i, c[q = 1.959 x 10~^ 



(7) 



and when q^ /ml = 0.5, one has. 



0.3125, Co = 1.1502, 



') "-2 



4 = 2.120 X 10"^ + 2.174 X lO^^z, 4 = -4.869 x 10"^ - 1.552 x 10 



-2 ■ 



4 = 1.420 X 10"^ + 5.174 X lO^^i, 4 = -5.729 x 10"^ - 1.552 x 10 



1-2' 



5 "-6 



4 = -8.340 X 10"^ - 9.938 x lO^^i, 4 = 3.839 x 10"^ 



-1.017 X 10"^ - 9.938 X 10"^^, 4o = 1.959 x 10"^ 



where we have taken as{mz) = 0.112, Q;em("^fe) = 1/132.2, nib = 5GeV, and rric = 
1.35GeV. 

3 CP violation in B^ -^ p^{uj)7T^ -^ tt^tt'tt^ 

3.1 Formalism 

The formahsm for CP violation in hadronic B meson decays is the following. Let A be 
the amplitude for the decay B^ — > 7r+7r^7r+, then one has 

A = {n+n-n+\H^\B+) + (7r+7r-7r+|i/^|5+), (9) 

with H^ and H^ being the Hamiltonians for the tree and penguin operators , respectively. 
We can define the relative magnitude and phases between these two contributions as 
follows, 

A = (7r+7r-7r+|//^|fi+)[l +re*V"^], (10) 

A = {n^TT-TT-\H^\B-)[l + re'^e-% (11) 

where 6 and (f) are strong and weak phases, respectively. The phase (j) arises from the 
appropriate combination of CKM matrix elements which is = a.Tg[{VtbV-i*j) / (VubV*^)]. 



As a result, sin0 is equal to sin a with a defined in the standard way |]Tl|]. The parameter 
r is the absolute value of the ratio of tree and penguin amplitudes: 



{p%u)7i+\H''\B+) 



{p^{uj)7T+\m\B+) 

The CP violating asymmetry, a, can be written as: 



(12) 



_|yl|2-|A|2_ -2r sin (5 sin (/) 
"" |A|2 + |A|2 " l + 2rcos(5cos0 + r2' ^^^^ 



It can be seen explicitly from Eq.(13) that both weak and strong phase differences are 
needed to produce CP violation. In order to obtain a large signal for direct CP violation, 
we need some mechanism to make both sin 5 and r large. We stress that p — uj mixing 
has the dual advantages that the strong phase difference is large (passing through 90° at 
the u resonance) and well known |Q, Q. With this mechanism, to first order in isospin 
violation, we have the following results when the invariant mass of vr+vr^ is near the uj 



(7r-7r+7r+|/7^|i?+) = ^^II,.t. + ^tp, (14) 



^p'^UJ ^p 



(7r-7r+7r+|iJ^|5+) = -^U^^p^ + ^p,. (15) 



°p°UJ °r 



Here tv{V = p or a;) is the tree amplitude and pv the penguin amplitude for producing 
a vector meson, V, Qp is the coupling for p° -^ vr+vr", Up^^ is the effective p — u mixing 
amplitude, and sy is from the inverse propagator of the vector meson V, 



sv = s -rriy + imvTv, (16) 



with -^i being the invariant mass of the vr+vr pair. 



We stress that the direct couphng uj -^ tt^tt is effectively absorbed into Up^j ^2\ 



leading to the explicit s dependence of Ilpuj. Making the expansion I[p^{s) = I[p^{m^) + 

(s — m^)n' (m^), the p — u mixing parameters were determined in the fit of Gardner and 

c>™ fr ('.^2\ _ _onn -u 'ir\r\T\/\a\i''^ 



O'Connell |T3l: ^e Ilp^{ml) = -3500 ± SOOMeV^ Sm Ilp^{ml) = -300 ± SOOMeV 
and il'pt^iml,) = 0.03 ± 0.04. In practice, the effect of the derivative term is negligible. 
From Eqs. (10,14, 15) one has. 






Defining 






where 6a, 6j3 and 5q are strong phases, one finds the following expression from Eq.(18) 



.^ _ ./.i5 n,. + /^e^^/'s. 



re'' = rV^^ ^ S J . (19) 



-'UJ \^ ^^pu)^ 



It will be shown that in the factorization approach, we have ae^^" = 1 in our case. Letting 



I3e'^p = b + ci, r'e'^" = d + ei, (20) 



and using Eq.(20), we obtain the following result when y/s ~ m^ 



,s C + Di 



re'^ = = ■ = , (21) 

(s - m^ + gfJe Up^y + {Qm Up^ + m^T^y 



where 



C={s-ml + ^e Up^)S^d[^e flp^ + b{s - ml) - cm^T^] 
-e[Qm f[p^ + bm^r^ + c{s - ml)] | 
+ (Sm Up^ + m^T^)le[^e lip^ + b{s - ml) - cm^T^] 
+d[Qm Ilp^ + bm^jT^ + c{s - ml)] L 

D = {s - ml + 'Re fi.puj)\ e[3?e 11^^ + d{s - ml) - cm^jT^] 
+d[^m Up^ + bm^T^ + c{s - ml)] \ 
-(3m fip^ + m^r^)ld[^e fip^ + b{s - ml) - cm^T^] 
-e[QmIlp^ + bm^T^ + c{s-ml)]y (22) 

Pe^^i^ and r'e'^^'' will be calculated later. Then, from Eq.(22) we can obtain rsin(5 and 
rcos6. In order to get the CP violating asymmetry, a, in Eq.(13), sin0 and cos0 
are needed, where is determined by the CKM matrix elements. In the Wolfenstein 



parametrization |M, one has 



sin0= , ' (23) 

cos0= , P(l-P)-^^ . (24) 



3.2 Calculational Details 

With the Hamiltonian given in Eq.(l), we are ready to evaluate the matrix elements for 
B^ —>■ p^ [00)71^. In the factorization approximation, either the p^{uj) or the tt^ is generated 
by one current which has the appropriate quantum numbers in the Hamiltonian. For this 
decay process, two kinds of matrix element products are involved after factorization; 
schematically (i.e. omitting Dirac matrices and color labels) (p°(u;)|(-un)|0)(7r+|(J6)|i?+) 
and {'n''^\{du)\0){p^{uj)\{ub)\B^). We will calculate them in some phenomenological quark 
models. 

The matrix elements ioi B ^ X and B -^ X* (where X and X* denote pseudoscalar and 
vector mesons, respectively) can be decomposed as [|T^, 

{X\MB) = L+px - ^^~r^ ^) F^iP) + "^^^^KF.iP), (25) 



TUB + nix* 

'' ■ ^ iPB + Px*),A2{e) - ^2mx* ■ k,A,{k')} 



niB + mx* k, 

+i^-^2mx* ■ k^Aoik^), (26) 
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where J^ is the weak current (J^ = g7^(l — 75)6 with q = u,d), k = ps — Px{x*) and 
e^ is the polarization vector of X*. The form factors included in our calculations satisfy: 
Fi(0) = Fo(0), ^3(0) = Ao(0) and Asik') = vi^i^A,{k^) - Vl^l^A^i^). Using the 
decomposition in Eqs.(25, 26), one has, 



tp = mslPp 



[c'l + ^J^2)fM^l) + (4 + ^c[)UA,{ml) 



(27) 



where fp and /^ are the decay constants of the p and tt, respectively, and Pp is the three 
momentum of the p. In the same way, we find t^i = tp, so that 



ae 



1. 



(2^ 



After calculating the penguin operator contributions, one has. 



/3e'' 



mslPp 



c'4 + Tr4)h/p^iK) + /.^oK)] 



N, 



+o[(4 + Tr4) + (4 + Tr4o)]/p^i(^ 



N, 



Nr, 



1 1 

-[(4 + ^4) + (4 + ^4)] 



2mlf^A^{ml 



{ttIu + md)imb + mu] 



+ (4o + ]^4)[^/p^iK) + UMml)]], 



r'e'^" 



Pu 



{c'l + ^/2)fpF,{ml) + (4 + ^/,)UAo{ml) 



VtbV*, 



VubV*d 



(29) 



where 



Puj = mB\Pp\{2 



(4 + Tr4 



Nr 



^^^ + iV^^« 



fpFiiml 
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1 

+ 2 
-2 



(4 + —4) + (4 + — c'lo) 



/pi^iK) 



m2/^Ao(m2) 



(m„ + md)(m(, + m„) 



+ (c^ + 



-4 



/,Ao(m2) + /^Fi(m' 



+ (cin + T^c' 



AT, 



^AoK)--/,Fi(m; 



and 



Vt,v:, 



VubV*a 






^~ 2 . 



sin7 



sin/? 



(30) 



3.3 Numerical Results 



In our numerical calculations we have several parameters: g^, N^. and the CKM matrix 
elements in the Wolfenstein parametrization. As mentioned in Section 2, the value of 
q^ is conventionally chosen to be in the range 0.3 < q^ jm^^ < 0.5. The CKM matrix, 
which should be determined from experimental data, has the following form in term of 



the Wolfenstein parameters. A, A, p, rj [M: 



V 



l-iA2 A AX^{p-i7]) 

-X 1 - iA2 AX' 

i3n ^ „•^^ A\2 I 



AX-^{1- p-i7]) -AX' 



(31) 



where O(A^) corrections are neglected. We use A = 0.2205, A = 0.815 and the range for 



p and 1] as the following [16, 17 



0.09 < p< 0.254, 0.323 < t] < 0.442. 



(32) 
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The form factors Fi(m^) and Ao{m1) depend on the inner structure of the hadrons. Under 
the nearest pole dominance assumption, the fc^ dependence of the form factors is: 



for model 1(2) pi, p 



F,{e 



h 



1- 



k2 ' 



Ao{k' 



hAo 



1-^' 



'Ao 



(33) 



where h = 0.330(0.625), /iaq = 0.28(0.34), mi = 5.32GeV, niAo = 5.27GeV, 



for model 3(4) p|, |T8|, |T|: 



F,{e 



h, 



k2 



2' 



Ao{k' 



*Ao 



1-:^ 



2' 



(34) 



where hi = 0.330(0.625), /i^o = 0.28(0.34), mi = 5.32GeV, m^o = 5.27GeV, 



for model 5 |20, 21|: 



F,ik' 



hi 



1 — ai 



fc2 



bi 



k'^ 



2' 



Ao(A;2 



/i 



Ao 



1 - aoK + bo{K 



2' 



(35) 



where hi = 0.305, /iaq = 0.372, ai = 0.266, 6i = -0.752, ao = 1.4, bo = 0.437. 
The decay constants used in our calculations are: fp = fu = 221MeV and /^r = 130.7MeV. 
In the numerical calculations, it is found that for a fixed A''^, there is a maximum value, 
o>maxi for the CP violating parameter, a, when the invariant mass of the tt+tt" is in the 
vicinity of the uj resonance. The results are shown in Figs.l and 2, for k'^/ml = 0.3(0.5) 
and Nf. in the range 0.98(0.94) < N^ < 2.01(1.95) - for reasons which will be explained 
later (Section 4). We investigate five models with different form factors to study the 
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model dependence of a. It appears that this dependence is strong (Table 1). 

The maximum asymmetry parameter, amax, varies from — 24%(— 19%) to — 59%(— 48%) 
for Nc in the chosen range, k"^ /ml = 0.3(0.5) and the range of CKM matrix elements in- 
dicated earlier. If we look at the numerical results for the asymmetries (Table 1) for 
Ncmax = 2.01(1.95) and k'^/ml = 0.3(0.5), we obtain for models 1, 3, 5 an asymmetry, 
ttmax, around — 27.3%(— 21.6%) for the set {pmax,Vmax), and around — 44.3%(— 35.0%) for 
the set {pmin,Vmin)- Wc find a ratio equal to 1.62(1.62) between the asymmetries associ- 
ated with the upper and lower limits of (p, rj). The reason why the maximum asymmetry, 
o-max, can have large variation, comes from the h ^ d transition, where Vm and Vuh appear. 
These are functions of (p, r^) and contribute to the asymmetry (Eq.31) through the ratio 
between the uo penguin diagram and the p tree diagram. 

For models 2 and 4, one has a maximum asymmetry, a.maxi around — 37%(— 28%) for 
the set {pmax,Vmax) and around — 59%(— 46%) for the set {pmin,Vmin)- We find a ratio 
between the asymmetries equal to 1.59(1.64) in this case. The difference between these 
two sets of models comes from the magnitudes of the form factors, where Fi{k'^) is larger 
for models 2 and 4 than for models 1, 3 and 5. Now, if we look at the numerical results 
for the asymmetry for Ncmin = 0.98(0.94), we find, for models 1, 3, 5, k'^/ml = 0.3(0.5), 
and the set {pmax,Vmax), an asymmetry, amax, around — 31.3%(— 25.6%), and for the set 
{pmin,Vmin) wc find an asymmetry, ttmax, around — 50.3%(— 42.0%). In this case, one has 
a ratio equal to 1.61(1.64). Finally, for models 2 and 4, we get — 36%(— 29%) for the set 
{Pmax,Vmax) and — 57%(— 48%) for the set {pmin,Vmin) with a ratio equal to 1.58(1.65). 

These results show explicitly the dependence of the CP violating asymmetry on form 
factors, CKM matrix elements and the effective parameter N^.. For the CKM matrix 

14 



elements, it appears that if we take their upper hmit, we obtain a smaller asymmetry, a, 
and viceversa. The difference between k"^ /ml = 0.3(0.5) in our results comes from the 
renormalization of the matrix elements of the operators in the weak Hamiltonian. Finally, 
the dependence on A''^ comes from the fact that A''^ is related to hadronization effects, and 
consequently we cannot determine Nc exactly in our calculations. Therefore, we treat 
Nc as a free effective parameter. As regards the ratio between the asymmetries, we have 
found a ratio equal to 1.61(1.63). This is mainly determined by the ratio sin7/sin/3, and 
more precisely by rj. In Table 2, we show the values for the angles a, f3, 7. From all 
these numerical results, we can conclude that we need to determine the value of A'c and 
the hadronic decay form factors more precisely, if we want to use the asymmetry, a, to 
constrain the CKM matrix elements. 

In spite of the uncertainties just discussed, it is vital to realize that the effect of 
p — uj mixing in the B ^ pn decay is to remove any ambiguity concerning the strong 
phase, sin 5. As the internal top quark dominates the b —>■ d transition, the weak phase 



VtdV* 



, and 



in the rate asymmetry is proportional to sin a (= sin0), where a = arg 
knowing the sign of sin 6 enables us to determine that of sin a from a measurement of the 
asymmetry, a. We show in Fig. 3 that the sign of sin 5 is always positive in our range, 
0.98(0.94) < Nc < 2.01(1.95) for all the models studied. Indeed, at the tt+tt" invariant 
mass where the asymmetry parameter, a, reaches a maximum, the value of sin 5 is equal 
to one - provided p — u mixing is included - over the entire range of N^ and for all the 
form factors studied. So, we can remove, with the help of asymmetry, a, the uncertainty 
mod(7r) which appears in a from the usual indirect measurements [^ which yield sin 2a. 
By contrast, in the case where we do not take p — uj mixing into account, we find a small 
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value for sin 5. In Figs. 3 and 4 we plot the role oi p — u mixing in our calculations. We 
stress that, even though one has a large value of sin 6 around Nc = 1 with no p — uj mixing, 
one still has a very small value for r (Fig. 4), and hence the CP violating asymmetry, a, 
remains very small in that case. 

4 Branching ratios for B^ -^ p^n^ and B^ -^ p^Ti~ 
4.1 Formalism 



With the factorized decay amplitudes, we can compute the decay rates using by the 



following expression |T9| 



T{B -^ VP) 



\Pp\ 



inrriy 



A{B -^ VP) 



e-pB 



(36) 



where 



\Pp\ 



J["m?B ~ ("^1 + ^712)^] [m^ — {nil ~ ""^2)^] 



2mB 



(37) 



is the cm. momentum of the decay particles, mi{m2) is the mass of the vector (pseu- 
doscalar) V(P), and A{B — > VP) is the decay amplitude: 



G 



A{B ^VP) = ^Y: Vj'^a,{VP\0,\B). 
V ^ i=i ,10 



(38) 



16 



Here Vj'^ is CKM factor: 



V:[=\VutV:,\ for t = l,2 and Vf 



\Vt,V,:i\ for ^ = 3,---,10 



where the effective parameters are the following combinations 



1 , 



1 , 



0'2j — C2j + -TTC2J-I1 (^2j-l — %-l + 1^%' fo^ J — I5 ■ ■ ■ 5 5 



and {VP\Oi\B) is a matrix element which is evaluated in the factorization approach. 
In the Quark Model, the diagram coming from the B^ -^ p^n^ decay is the only one 
contribution. In our case, to be consistent, we should also take into account the p — uo 
mixing contribution when we calculate the branching ratio since we are working to the 
first order of isospin violation. Explicitly, we obtain for B^ -^ p^n^, 



BR{B+ -^ p%+) 



Gl\p, 



V^A^o (ai, aa) - V^A^o (ag, ■ ■ ■ , oio) 



+ 



327rrB 

T aT/„ „ N T/P aP 



VjAUa,,a,)-V:'A^{a,, 



< o-io) 



n 



pUJ 



{sp - ml) + im^T^ 



(39) 



where the tree and penguin amplitudes are: 



\/2Aj)(ai,a2) = aifpFi{ml) + a2UAo{ml), 
V2Apo{a3, ■ ■ ■ , aio) = 04 -fpFi{ml) + f^Ao{ml 



+ aw 



+ -{aT + ag)fpFi{ml) - 2{aQ + as) 
V2Al{ai,a2) = aifpFi{ml) + a2UAo{ml), 



-fpF,{ml) + UA,{ml 
{iTLu + md){mi, + rriu) 
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V^^<^(«3, ■ ■ ■ , Oio) 



2(03 + 05) + -{a7 + ag) 



fpFMD 



-2(a8 + oe) 



mlUAo{ml) 



{ruu + md){mb + mu) 



+04 



UAoiml) + f^F,{m 



+ a 



10 



UAoiml) - -f,F,{ml) 



where (p°|mm|0) = -j^fpUipep and {7r+\ud\0) = if^^p^. 



For 5° -^ p+TT we obtain, 



5i?(5° -^ p+TT" 



IGvrr^o 



If Aj+(a2) - Kf <+(«3, ■ ■ ■ , Oio) 



(40) 



where 



^^+(02) =a2fpFi{ml), 
^r+(«3,---,aio) = (a4 + aio)/pi^i("^p)- 



Moreover, we can calculate the ratio between these two branching ratios, in which the 
uncertainty caused by many systematic errors is removed. We define the ratio R as: 



R 



BR{B^ -^ p+TT-) 
BR{B+ ^ pH+y 



(41) 



and, without taking into account the penguin contribution, one has. 



R 



2Vb+ 



r_Bo 



02 fpFiiml) 



n 



piii 



[So -mi) + 



im^^^ 



(42) 
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4.2 Numerical Results 

The latest experimental data from the CLEO collaboration ^ are: 

BR{B+ -^ p°7r+) = (10.4+1^ ± 2.1) x 10-^ 
BR{B^ -^ p+TT") = (27.6+?i ± 4.2) x 10-^ 
i? = 2.65 ±1.9. 

We have calculated the branching ratios for B^ -^ p+vr^ and for B^ -^ p°7r+ for all models 
as a function of Nc. In Figs. 5 and 6, we show the results for models 1 and 2 in order to 
make the dependence on form factors explicit. 

The numerical results are very sensitive to uncertainties coming from the experimental 
data. For the branching ratio B^ -^ p+7r~ (Fig.5), we have a large range of values of Nc 
and the CKM matrix elements over which the theoretical results are consistent with the 
experimental data from CLEO. However, all models do not give the same result: models 
2 and 4 are very close to the experimental data for a large range of N^ whereas models 

1.3 and 5 are not. The reason is still the magnitude of the form factors. As a result, we 
have to exclude models 1,3 and 5 because their form factors are too small. 

If we consider numerical results for branching ratio B^ -^ p^vr"*" (Fig.6), it appears that 
all models are consistent with the experimental data for a large range of Nf.. The effect of 
p — uo mixing (included in our calculations) on the branching ratio B^ — > p^vr^ is around 
30%. Numerical results for models 1, 3, 5 and models 2, 4 are very close to each other. 
The difference between the two branching ratios can be explained by the fact that for 
the B^ -^ p^-K^ decay, the tree and penguin contributions are both proportional to only 
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one form factor, Fi{k'^). Thus, this branching ratio is very sensitive to the magnitude 
of this form factor (Fi(A;^) is related to hi = 0.330 or 0.625 in models (1,3) and (2,4) 
respectively). On the other hand, for the decay B^ -^ P^^r^, both Fi{k'^) and Ao[k'^) are 
included in the tree and penguin amplitudes, and this branching ratio is less sensitive to 
the magnitude of the form factors. 

If we look at the ratio R between these two branching ratios, BR{B^ -^ p^ii^) and 
BR{B^ -^ P^^ ) ^ shown in Fig. 7 - the results indicate that R is very sensitive to the 
magnitude of the form factors, and that there is a large difference between models 1,3,5 
and models 2 and 4. We investigated the ratio R for the limiting CKM matrix elements 
as a function of Nc, finding that R is consistent with the experimental data over the range 
Nc- 0.98(0.94) < Nc < 2.01(1.95), (The values outside (inside) brackets correspond to the 
choice q^ /ml = 0.3(0.5)). It should be noted that R, in particular, is not very sensitive 
to the CKM matrix elements. The small difference which does appear, comes from the 
penguin contributions (which may be neglected). If we just take into account the tree 
contributions in our calculations, R is clearly independent of the CKM matrix elements 
(Eq.42). 

From a comparison of the numerical results and the experimental data, we can extract 
a range of N^, within which all results are consistent. In Table 3, we have summarized 
the allowed range of Nc for B^ -^ P°7r+, B^ -^ p+vr^ and R, for models 1, 2, 3, 4 and 5 
according to various choices of the CKM matrix elements. To determine the best range of 
Nc, we have to find some intersection of the values of Nc for each model and for each set of 
CKM matrix elements, for which the theoretical and experimental results are consistent. 
This is possible and the results are shown in Table 4. In our study, it seems better to use 
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the range intersection {Nc}b+ H {Nc}r than {Nc}bo H {Nc}b+, for fixing the final interval 
Nc, since the experimental uncertainties are smaller in the former case, and since we 
are working to the first order of isospin violation {p — uj mixing). Finally, after excluding 
models 1,3 and 5, which are not consistent with all the experimental data, we are able to fix 
the upper and lower limit of the range of Nc, using the limiting values of the CKM matrix 
elements (Table 5). We find that N^ should be in the range 0.98(0.94) < N^ < 2.01(1.95) 
where Ncmin and Ncmax correspond to {pmin,Vniin) and {pmax,Vniax) respectively. 

5 Summary and discussion 

The first aim of the present work was to compare our theoretical results with the latest 
experimental data from the CLEO collaboration for the branching ratios B^ -^ p°7r+ 
and B^ -^ p+vr^. Our next aim was to study direct CP violation for the decay B^ -^ 
p^{(jj)-n^ -^ 7r+7r^7r+, with the inclusion of p — a; mixing. The advantage of p — tu mixing 
is that the strong phase difference is large and rapidly varying near the uo resonance. As 
a result the CP violating asymmetry, a, has a maximum, amax) when the invariant mass 
of the TT+TT" pair is in the vicinity of the u resonance and sin 5 = +1 at this point. 

In the calculation of CP violating asymmetry parameters, we need the Wilson coeffi- 
cients for the tree and penguin operators at the scale m^. We worked with the renormal- 
ization scheme independent Wilson coefficients. One of the major uncertainties is that 
the hadronic matrix elements for both tree and penguin operators involve nonperturbative 
QCD. We have worked in the factorization approximation, with Nc treated as an effective 
parameter. Although one must have some doubts about factorization, it has been pointed 
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out that it may be quite reliable in energetic weak decays |2^, ^ . 

We have explicitly shown that the CP violating asymmetry, a, is very sensitive to the 
CKM matrix elements and the magnitude of the form factors, and we have determined a 
range for the maximum asymmetry, ttmaxi as a function of the parameter Nc-, the limits 
of CKM matrix elements and the choice of k'^/ml = 0.3(0.5). From all the models 
investigated, we found that CP violating asymmetry, amax, varies from — 24%(— 19%) to 
— 59%(— 48%). We stressed that the ratio between the asymmetries associated with the 
limiting values of CKM matrix elements would be mainly determined by t]. Moreover, we 
also stressed that without p — u mixing, we cannot have a large CP violating asymmetry, 
a, since a is proportional to both sin6 and r. Even though sin (5 is large around Nc = 1, 
r is very small. As a result, we find a very small value for the CP violation in the decay 
B^ -^ p^TT^ (of the order of a few percent) without mixing. Once mixing is included, 
the sign of sin 5 is positive for N^ : 0.98(0.94) < N^ < 2.01(1.95). Indeed, at the tt+tt' 
invariant mass where the asymmetry, a, is maximum, sin 5 = +1, independent of the 
parameters used. Thus, by measuring a, we can erase the phase uncertainty mod(7r) in 
the determination of the CKM angle a which arises from the conventional determination 
of sin 2a. 

The theoretical results for the branching ratios for B^ -^ p^ir^ and 5° -^ P^t^^ , were 
compared with the experimental data from the CLEO collaboration 0. These calculations 
show that it is possible to have theoretical results consistent with the experimental data 
without needing to invoke contributions from other resonances P^, |23|. These data helped 



us to constrain the magnitude of the various form factors needed in the theoretical calcu- 
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lations of B decays^. We determined a range of value of Nc, 0.98(0.94) < Nc < 2.01(1.95), 
inside of which the experimental data and the theoretical calculations are consistent for 
models 2 and 4. 

We will need more accurate data in the future to further decrease the uncertainties in 
the calculation. If we can use both the CP violating asymmetry and the branching ratios, 
with smaller uncertainties, we expect to be able to determine the CKM matrix elements 
more precisely. At the very least, it appears that one will be able to unambiguouly de- 
termine the sign of sin a and hence, remove the well known discrete uncertainties in a 
associated with the fact that indirect CP violation determines only sin 2a. We expect 
that our predictions should provide useful guidance for future investigations and urge 
our experimental collegues to plan seriously to measure the rather dramatic direct CP 
violation predicted here. 
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Figure Captions 

Fig.l Asymmetry, a, for k'^/ml = 0.3, Nc = 0.98(2.01) and limiting values of the CKM 
matrix elements for model 1: solid line(dot line) for N^ = 0.98 and max(min) CKM ma- 
trix elements. Dashed line(dot dashed line) for Nc = 2.01 and max(min) CKM matrix 
elements. 

Fig. 2 Asymmetry, a, for k'^/ml = 0.5, N^ = 0.94(1.95) and limiting values of the CKM 
matrix elements for model 1: solid line(dot line) for A*",, = 0.94 and max(min) CKM ma- 
trix elements. Dashed line(dot dashed line) for Nc = 1.95 and max(min) CKM matrix 
elements. 

Fig. 3 Determination of the strong phase difference, sin 5, for k'^/m\ = 0.3(0.5) and for 
model 1. Solid line(dot hne) for Up^^ = (—3500; —300) (i.e. with p — u mixing). Dot 
dashed line(dot dot dashed line) for IIp^ = (0; 0), (i.e. with no p — u mixing). 

Fig. 4 Evolution of the ratio of penguin to tree amplitudes, r, for k'^/m'^ = 0.3(0.5), for lim- 
iting values of the CKM matrix elements (p, rj) max(min), for Ilp^ = (—3500; — 300)(0, 0), 
(i.e. with(without) p — u mixing) and for model 1. Figure 4a (left): for Up^^ = (0;0), 
solid line(dot line) for k'^/m\ = 0.3 and {p,i]) max(min). Dot dashed line(dot dot dashed 
line) for k'^/m?^ = 0.5 and {p,ri) max(min). Figure 4b (right): same caption but for 
n^^ = (-3500; -300). 
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Fig. 5 Branching ratio for B^ -^ p+vr^ for models 1(2), k'^/m\ = 0.3 and limiting val- 
ues of the CKM matrix elements. Solid line(dot line) for model 1 and max(min) CKM 
matrix elements. Dot dashed line (dot dot dashed line) for model 2 and max(min) CKM 
matrix elements. 

Fig. 6 Branching ratio for B^ -^ p°7r+ for models 1(2), k'^/m\ = 0.3 and limiting val- 
ues of the CKM matrix elements. Solid line(dot line) for model 1 and max(min) CKM 
matrix elements. Dot dashed line (dot dot dashed line) for model 2 and max(min) CKM 
matrix elements. 

Fig. 7 Calculation of the ratio of the two pvr branching ratios versus Nc for models 1(2) 
and for limiting values of the CKM matrix elements: solid line(dot line) for model 1 with 
max(min) CKM matrix elements. Dot dashed line(dot dot dashed line) for model 2 with 
maxfmin) CKM matrix elements. 
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Tables 

Table 1 Maximum CP violating asymmetry, amax{%), for B^ -^ tt^tt^tt^, for all models, 
limiting values of the CKM matrix elements (upper and lower limit), and for k"^ /ml = 
0.3(0.5). 

Table 2 Values of the CKM unitarity triangle for limiting values of the CKM matrix 
elements. 

Table 3 Summary of the range of values of Nc which is determined from the experi- 
mental data for various models and input parameters. 

Table 4 Determination of the intersection of the values of Nc which are consistent with 
various subsets of the data for all models and all sets of CKM matrix elements. 

Table 5 Best range of Nc determined from Table 4 for k'^/ml = 0.3(0.5). 
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Figure 1: Asymmetry, a, for k'^ /ml = 0.3, N^ = 0.98(2.01) and limiting values of the 
CKM matrix elements for model 1: solid line(dot line) for Nc = 0.98 and max(min) CKM 
matrix elements. Dashed line(dot dashed line) for N^. = 2.01 and max(min) CKM matrix 
elements. 
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Figure 2: Asymmetry, a, for k"^ /ml = 0.5, Nc = 0.94(1.95) and limiting values of the 
CKM matrix elements for model 1: solid line(dot line) for Nc = 0.94 and max(min) CKM 
matrix elements. Dashed line(dot dashed line) for Nc = 1.95 and max(min) CKM matrix 
elements. „„ 




Figure 3: Determination of the strong phase difference, sin 5, for k'^ /m\ = 0.3(0.5) and for 
model 1. The sohd(dotted) hne at sin 5 = +1 corresponds the case Ilp^ = (—3500; —300), 
where p — u mixing is included. The dot dashed(dot dot dashed) line corresponds to 
Hpw = (0; 0), where p — co mixing is not included. 
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Figure 4: Evolution of the ratio of penguin to tree amplitudes, r, for k'^ /m\ = 
0.3(0.5), for limiting values of the CKM matrix elements (p, //) max(min), for Up^^ = 
(—3500; — 300)(0, 0), (i.e. with(without) p — u mixing) and for model 1. Figure 4a (left): 
for ITp^ = (0;0), solid line(dot line) for k'^/m% = 0.3 and (p, 77) max(min). Dot dashed 
line(dot dot dashed line) for k'^/m\ = 0.5 and (p, //) max(min). Figure 4b (right): same 
caption but for Up^ = (—3500; —300). 39 
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Figure 5: Branching ratio for B^ -^ p+vr^ for models 1(2), k'^/m'^ = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dot line) for model 1 and max(min) CKM 
matrix elements. Dot dashed line (dot dot dashed line) for model 2 and max(min) CKM 
matrix elements. 
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Figure 6: Branching ratio for B~^ —* p^ir^ for models 1(2), k'^/m% = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dot line) for model 1 and max(min) CKM 
matrix elements. Dot dashed line (dot dot dashed line) for model 2 and max(min) CKM 
matrix elements. 
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Figure 7: Calculation of the ratio of two pvr branching ratios versus N^. for models 1(2) 
and for limiting values of the CKM matrix elements: solid line(dot line) for model 1 with 
max(min) CKM matrix elements. Dot dashed line(dot dot dashed line) for model 2 with 
maxfmin) CKM matrix elements. 
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iV,„,„ = 0.98(0.94) 


iV,„,, = 2.01(1.95) 


model 1 


pmax 1 Vmax 
Pmin 1 Vmin 


-33(-27) 
-52 (-43) 


-29(-23) 

-47(-37) 


model 2 


Pmax 1 Vmax 
Pmin 1 Vmin 


-36(-29) 

-57(-48) 


-37(-28) 
-59(-46) 


model 3 


Pmax 1 Vmax 
pmin 1 Vmin 


-32(-26) 
-51 (-43) 


-29(-23) 
-47(-37) 


model 4 


Pmax 5 Vmax 
Pmin ; Vmin 


-36(-29) 

-57(-48) 


-37(-28) 
-59(-46) 


model 5 


Pmax ) Vmax 
Pmin ; Vmin 


-29 (-24) 
-48 (-40) 


-24(-19) 
-39(-31) 



Table 1: Maximum CP violating asymmetry amax{%) for B~^ -^ vr+vr vr"*", for all models, 
limiting values of the CKM matrix elements (upper and lower limit), and for k'^/ml = 
0.3(0.5). 





ip,V)min 


(,P) V jniax 


a 


86°02 


89°23 


P 


19^50 


30°64 


7 


74043 


60°11 



Table 2: Values of the CKM unitarity triangle for limiting values of the CKM matrix 
elements. 
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5+ 5° R 


model 1 


Pmax^riraax 0.76; 1.69(0.73;1.62) 5.50; ^^ ( - ; - ) 0.92;2.57(0.90;2.52) 


Pmin^rimin 0.52; 1.04(0.49;0.98) 


-(- 


- ) 0.97;2.88(0.94;2.76) 


Pma.,Vm^n 0.61; 1.25(0.59;1.20) 


-(- 


- ) 0.92;2.58(0.91;2.54) 


Prmn,V^ax 0.69; 1.46(0.66;1.39) 


-(- 


- ) 0.95;2.75(0.90;2.66) 


model 2 


pmax,Vmax 1.44;3.06(1.40;2.95) 0.54;1.33(0.54;1.38) 0.86;1.89(0.84;1.86) 


pmin,Vmin 1.00;2.01 (0.96;1.90) 1.10; ^ (1.15; ^*) 0.92;2.09(0.89;2.01) 


pmax,Vmin 1.15;2.32(1. 12;2.22) 0.70; ^ (0.72; ^* ) 0.87;1.89(0.85;1.86) 


Prmn,Vmax 1.32;2.78(1.25;2.60) 0.63;2.77(0.62;3.12) 0.90;2.00(0.84;1.94) 


model 3 


Pmax,Vn.ax 0.74; 1.65(0.72;1.60) 


-(- 


- ) 0.92;2.65(0.92;2.60) 


pmin,Vmin 0.51; 1.02(0.49;0.98) 


-(- 


- ) 0.97;2.95(0.94;2.85) 


Pmax,Vm^n 0.60; 1.22(0.57;1.19) 


-(- 


- ) 0.93;2.66(0.92;2.61) 


pmin,Vmax 0.67;1.43(0.65;1.37) 


-(- 


- ) 0.92;2.79(0.92;2.71) 


model 4 


pmax,Vniax 1.41;3.04(1.36;2.92) 0.56;1.44(0.57;1.52) 0.86;1.91(0.85;1.87) 


pmin,Vmin 0.98; 1.96(0.94;1.87) 1.16; ^ (1.23; ^^ ) 0.90;2.10(0.89;2.03) 


Pmax,Vrmn 1.14;2.29(1. 10;2.21) 0.72; ^ (0.74; ^^ ) 0.86;1.92(0.85;1.88) 


pmin,Vmax 1.30;2.74(1.24;2.59) 0.64;3.49(0.66;4.03) 0.89;2.01(0.86;1.95) 


model 5 


Pmax 1 Vmax U . ( 


2.18(0.73;2.10) 


-(- 


- ) 1.03; Tb^ (1.02; i.i. ) 


PminyVmin U.OU 


1.08(0.47;1.03) 


-(- 


- ) 1.09; ^ (1.06; ^^ ) 


PmaxyVmin U.Oo 


1.38(0.55;1.34) 


-(- 


- ) 1.03; ^ (1.02; ^* ) 


Pmini Vmax U.OO 


1.71(0.64;1.62) 


-(- 


- ) 1.04; ^ (1.04; *^ ) 



Table 3: Summary of the range of values of N^ which is determined from the experimental 
data for various models and input parameters (numbers outside(inside) brackets are for 
k^/ml = 0.3(0.5)). The notation: 

(number ; number) means that there is a upper and lower limit for A^^- (number ; -k-k ) 
means that there is no upper limit for Nc in the range Nc [0;10]. ( - ; - ) means that 
there is no range of N^ which is consistent with experimental data. 
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{n,}b+ n {n,}bo {ivj^+ n {n,}r {n,}bo n {n,}r 


model 1 


pmax,Vmax -{-) 0.92; 1.69(0.90;1.62) -(-) 


Pmin,Vn.in -(-) 0.97; 1.04(0.94;0.98) -(-) 


Pmax,Vmin -(-) 0.92; 1.25(0.91;1.20) -(-) 


pmin,Vmax -{-) 0.95; 1.46(0.90;1.39) -(-) 


model 2 


Pmax,Vmax -(-) 1.44; 1.89(1.40;1.86) 0.86;1.33(0.84;1.38) 


pmin,Vmin 1.10;2.01 (1. 15;1.90) 1.00;2.01(0.96;1.90) 1.10;2.09(1.15;2.01) 


Pmax,Vmin 1.15;2.32(1. 12;2.22) 1.15;1.89(1.12;1.86) 0.87;1.89(0.85;1.86) 


Prmn,Vmax 1.32;2.78(1.25;2.60) 1.32;2.00(1.25;1.94) 0.90;2.00(0.84;1.94) 


model 3 


pmax,Vrr.ax -{-) 0.92;1.65(0.92;1.60) -(-) 


pmin,Vmin -{-) 0.97; 1.02(0.94;0.98) - ( - ) 


Pmax,Vmin -(-) 0.93; 1.22(0.92;!. 19) -(-) 


Pmin,Vmax -(-) 0.92; 1.43(0.92;1.37) -(-) 


model 4 


pmax,Vmax 1.41; 1.44(1. 36;1.52) 1.41;1.91(1.36;1.87) 0.86;1.44(0.85;1.52) 


pmin,Vmin 1.16; 1.96(1. 23;1.87) 0.98;1.96(0.94;1.87) 1.16;2.10(1.23;2.03) 


pmax,Vmin 1.14;2.29(1. 10;2.21) 1.14;1.92(1.10;1.88) 0.86;1.92(0.85;1.88) 


pmin,Vmax 1.30;2.74(1.24;2.59) 1.30;2.01(1.24;1.95) 0.89;2.01(0.86;1.95) 


model 5 


Pn^ax,Vmax -(-) 1.03;2. 18(1.02;2.10) -(-) 


Pmin 5 Vmin ^ \ ^ ) ^ (, ^ ) ^ (, ^ ) 


Pmax,Vrmn -{-) 1.03; 1.38(1.02;1.34) -(-) 


Prmn,V^ax -{-) 1.04; 1.71 (1.04;1.62) -(-) 



Table 4: Determination of the intersection of the values of N^. which are consistent with 
various subsets of the data for all models and all sets of CKM matrix elements (numbers 
outside(inside) brackets are for k'^/ml = 0.3(0.5)). The notation: 
- ( - ) means that no common range of N^ can be extracted from the data. 
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{Nc} with mixing {Nc} without mixing 

model 2 1.00;2.01(0.96;1.94) 0.85;1.74(0.85;1.74) 

model 4 0.98;2.01(0.94;1.95) 0.84;1.76(0.84;1.75) 

maximum range 0.98;2.01(0.94;1.95) 0.84;1.76(0.84;1.75) 

minimum range 1.00;2.01(0.96;1.94) 0.85;1.74(0.85;1.74) 



Table 5: Best range of Nc determined from Table 4 for k'^/ml = 0.3(0.5). One takes the 
maximum interval of Nc, from Table 4, for each model (2,4). To determine the maxi- 
mum(minimum) range, one consideres all models (2,4) and the largest (smallest) range of 
Nc- In comparison, we show the range of Nc determined without p — uj mixing. 
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